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Abstract 

I sketch the architecture of O 'Mega, a new optimizing compiler for 
tree amplitudes in quantum field theory. O'Mega generates the most 
efficient code currently available for scattering amplitudes for many 
polarized particles in the standard model. A complete infrastructure 
for physics beyond the standard model is provided. 



1 Introduction 

Current and planned experiments in high energy physics can probe processes 
with many tagged — potentially polarized — particles in the final state. The 
combinatorial explosion of the number of Feynman diagrams contributing 
to scattering amplitudes for many external particles calls for the develop- 
ment of more compact representations that translate well to efficient and 
reliable numerical code. In gauge theories, strong numerical cancellations in 
a redundant representation built from necessarily gauge dependent Feynman 
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diagrams lead to a loss of numerical precision, stressing further the need for 
eliminating redundancies. 

Due to the large number of processes that have to be studied in order 
to unleash the potential of modern experiments, the construction of these 
representations must be possible algorithmically on a computer and should 
not require human ingenuity for each new application. 

O'Mega lU is a compiler for tree- level scattering amplitudes that satisfies 
these requirements. O'Mega is independent of the target language and can 
support code in any programming language for which a simple output module 
has been written. To support a physics model, O'Mega requires as input only 
the Feynman rules and the relations among coupling constants. 

Similar to earlier numerical approaches 0, [|, O'Mega reduces the growth 
in calculational effort from a factorial of the number of particles to an ex- 
ponential. The symbolic nature of O'Mega, however, increases its flexibility. 
Indeed, O'Mega can emulate both 0, Q and produces code that is empirically 
at least twice as fast. 



2 lPOWs And Keystones 



One Particle Off-shell Wave functions (lPOWs) are obtained from Greens- 
functions by applying the LSZ reduction formula to all but one line: 
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The 1POW W^ q \x) = ((f)(q),c/)(q');out\&(x)\(j)(p);m) in lowest order of 3 - 
theory, is given — for illustration — by 
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At tree-level, the set of all lPOWs for a given set of external momenta can 
be constructed recursively ||] 






k+l=n 
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where the sum extends over all partitions of the set of n momenta. For all 
quantum field theories, there are — well defined, but not unique — sets of Key- 
stones K |1J such that the sum of tree Feynman diagrams can be expressed 
as a sparse sum of products of lPOWs without double counting. In a theory 
with only cubic couplings this is expressed as 

F{n) P(n) 

T =J2 D ^= E K lhf m (Pk,PhPm)W fk ( Pk )W ft ( Pl )W fm ( Pm ) . (4) 
i=l k,l,m=l 

The non-trivial problem of avoiding the double counting of diagrams like (the 
circle denotes the keystone) 





has been solved for general theories with vertices of arbitrary degrees [|TJ. 

The number of distinct momenta that can be formed from n external 
momenta is P(n) = 2 n ~~ 1 — 1. Therefore, the number of tree lPOWs grows 
exponentially with the number of external particles and not with a factorial, 
as the number of Feynman diagrams F{n) = (2n — 5)!! = (2n — 5) ■ ... 5 ■ 
3-1. The equations sketched in Eqs. ([|) and ([|) for cubic couplings can be 
generalized to vertices of any order Jl|]. 

Even for vector particles and to all orders in renormalized perturbation 
theory, the lPOWs are 'almost' physical objects and satisfy simple Ward 
identities in unbroken gauge theories 

(outl^x) lin* =0 (5) 
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and well as in spontaneously gauge theories 




Qx MIW^Namp. = iwm W (out | (M^) I m> amp. ( 6 ) 

in Z?£-gauge. The code for matrix elements can optionally be instrumented 
by O'Mega with numerical checks of these Ward identities for intermediate 
lines. 
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3 Directed Acyclical Graphs 

The algebraic expression for the tree-level scattering amplitude in terms of 
Feynman diagrams is itself a tree. The much slower growth of the set of 
lPOWs compared to the set of Feynman diagrams shows that this repre- 
sentation is extremely redundant. In this case, Directed Acyclical Graphs 
(DAGs) provide a more efficient representation, as illustrated by a trivial 
example 



where one multiplication is saved. The replacement of expression trees by 
equivalent DAGs is part of the repertoire of optimizing compilers, known 
as common subexpression elimination. Unfortunately, this approach fails 
for typical expressions appearing in quantum field theory, because of the 
combinatorial growth of space and time required to find an almost optimal 
factorization. 

However, the recursive definition in Eq. (|3|) allows to construct the DAG of 
the lPOWs in Eq. (|J) directly U, without having to construct and factorize 
the Feynman diagrams explicitely. 

As mentioned above, there is more than one consistent prescription for 
constructing the set of keystones . The symbolic expressions constructed 
by O'Mega contain the symbolic equivalents of the numerical expressions 
computed by (maximally symmetric keystones) and H (maximally asym- 
metric keystones) as special cases. 

4 Algorithm 

By virtue of their recursive construction in Eqs. (^), tree- level lPOWs form 
a DAG and the problem is to find the smallest DAG that corresponds to a 
given tree, (i.e. an given sum of Feynman diagrams). O'Mega's algorithm 
proceeds in four steps 

Grow: starting from the external particles, build the tower of all lPOWs 
up to a given height (the height is less than the number of exter- 
nal lines for asymmetrical keystones and less than half of that for 
symmetrical keystones) and translate it to the equivalent DAG D. 
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Select: from D, determine all possible flavored keystones for the pro- 
cess under consideration and the lPOWs appearing in them. 

Harvest: construct a sub-DAG D* C D consisting only of nodes that 
contribute to the lPOWs appearing in the flavored keystones. 

Calculate: multiply the lPOWs as specified by the keystones and sum 
the keystones. 

By construction, the resulting expression contains no more redundancies and 
can be translated to a numerical expression. In general, asymmetrical key- 
stones create an expression that is smaller by a few percent than the result 
from symmetrical keystones, but it is not yet clear which approach produces 
the numerically more robust results. 

5 Implementation 

The O'Mega compiler is implemented in O'Caml ||, a functional program- 
ming language of the ML family with a very efficient, portable and freely 
available implementation, that can be bootstrapped on all modern comput- 
ers in a few minutes. 

The powerful module system of O'Caml allows an efficient and concise 
implementation of the DAGs for a specific physics model as a functor appli- 
cation jlj. This functor maps from the category of trees to the category of 
DAGs and is applied to the set of trees defined by the Feynman rules of any 
model under consideration. 

The implementation is concise and efficient simultaneously by exploiting 
the virtues of persistent data structures 0. Typically, the resources con- 
sumed by O'Mega are only a small fraction of the resources required by the 
compiler for the target language. 

The module system of O'Caml has been used to make the combinato- 
rial core of O'Mega demonstrably independent from the specifics of both 
the physics model and the target language [fT), as shown in Figure [I]. A 
Fortran90/95 backend has been realized first, backends for C++ and Java 
will follow. The complete electroweak standard model has been implemented 
(the treatment of interfering color amplitudes is still incomplete). Majorana 
fermions, required by supersymmetric field theories, are available (using M) 
and the MSSM is in preparation. 
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Figure 1: Module dependencies in O'Mega. The diamond shaped nodes 
denote abstract signatures defining functor domains and co-domains. The 
rectangular boxes denote modules and functors, while oval boxes stand for 
example applications. 
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As mentioned above, the compilers for the target programming language 
are the slowest step in the generation of executable code. On the other 
hand, the execution speed of the code is limited by non-trivial vertex eval- 
uations for vectors and spinors, which need O(10) complex multiplications. 
Therefore, an O'Mega Virtual Machine can challenge native code and avoid 
compilations. 



6 Applications 

The code generated by the Fortran90/95 backend is the most efficient code 
available for polarized scattering amplitudes for many particles. The results 
have been compared with MADGRAPH for many standard model pro- 
cesses and numerical agreement at the level of 10 -11 has been found with 
double precision floating point arithmetic. O'Mega generated amplitudes are 
used in the omnipurpose event generator generator WHIZARD ||. The first 
complete experimental study of vector boson scattering in six fermion pro- 
duction for linear collider physics [|l(J has been facilitated by O'Mega and 
WHIZARD. 
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